b = np.insert(a, insert_index, values=a[:,2], axis=1)
import numpy as np N = 10 a = np.random.rand(N,N) b = np.zeros((N,N+1)) b[:,:-1] = a